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We present the final results of the CP-PACS calculation of the light hadron spectrum and quark masses with 
two flavors of dynamical quarks. Simulations are made with a renormalization-group improved gauge action and 
a mean-field improved clover quark action for sea quark masses corresponding to mps/mv ~ 0.8-0.6 and the 
lattice spacing a = 0.22-0.11 fm. For the meson spectrum in the continuum limit a clearly improved agreement 
with experiment is observed compared to the quenched case, demonstrating the importance of sea quark effects. 
For light quark masses we obtain m^(2GeF) = 3.44t ;^ MeV and mY I {2GeV) = 88±g MeV (if -input) and 
m,g IS (2GeV) = 90!tii MeV (<?i-hiput), which are reduced by about 25% compared to the values in quenched QCD. 



1. INTRODUCTION 

Recent progress in lattice calculations of the 
quenched light hadron spectrum [jl| has demon- 
strated clearly that effects of dynamical sea 
quarks have to be included for a precision com- 
parison with the experimental spectrum. As a 
first step towards the real world we have pursued 
a systematic calculation in two-flavor QCD over 
the last 3 years §||. In this article we present a 
summary of the final results on the light hadron 
spectrum and quark masses from this calculation. 

In order to deal with the increased demand in 
computer time for full QCD we have used an RG- 
improved gauge action and a tadpole-improved 
SW clover quark action, which allows us to simu- 
late at coarse lattices in the range a = 0.22-0.11 
fm m. Since Lattice'99 [|| we have doubled 
the statistics at the smallest lattice spacing on 
a 24 3 x 48 lattice, with which we have reached 
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our planned target of three lattice spacings, each 
with a similar amount of statistics at each sim- 
ulated sea quark mass. The final parameters of 
our simulations are summarized in Table |l|. 

We have also made a set of simulations in 
quenched QCD using the same improved actions 
in the same range of lattice spacings and with 
a similar range of quark masses. This allows 
us to compare in detail spectrum results in full 
and quenched QCD and to make an independent 
confirmation of the findings in quenched QCD in 
Ref. 0] . Parameters of quenched simulations are 
summarized in Table ^. 

2. SEA QUARK EFFECTS IN THE 
LIGHT HADRON SPECTRUM 

A novel feature in the spectrum analysis with 
two-flavor QCD is the dependence of hadron 
masses on both the sea and the valence quark 
mass. We take this into account by a combined 
quadratic ansatz. For vector mesons and decuplet 
baryons the formula has the form, 

m>v,D = A + B s /j, sea + B v iJ, va i 

+Cs/J-lea + Cv^lal + CsvUseaHval , (1) 



2 



Table 1 

Parameters and final statistics of two-flavor QCD simulations. 





L 6 x T 


csw 


a[fm] 


-La[fm] 


mps/my for sea quarks 


: # HMC trajectories 




1.80 
1.95 
2.10 


12 3 x24 
16 3 x32 
24 3 x48 


1.60 
1.53 
1.47 


0.215(2) 
0.155(2) 
0.108(1) 


2.58(3) 
2.48(3) 
2.58(3) 


0.807(1):6250 0.753(1):5000 
0.804(1):7000 0.752(1):7000 
0.806(1):4000 0.755(2):4000 


0.694(2):7000 0.547(4) 
0.690(1):7000 0.582(3) 
0.691(3):4000 0.576(3) 


5250 
5000 
4000 



where \i va i represents the mass squared m 2 PS of 
a pseudo scalar meson made of a valence quark- 
antiquark pair and /x se a stands for a PS meson 
made of sea quarks. Since octet baryon masses 
are not functions of the average valence quark 
mass we use a more complicated formula with 
12 parameters. For quenched results with im- 
proved actions terms in /i sea are dropped and 
we find that vector meson and decuplet baryon 
masses are well described by a linear ansatz in 
Hvai- Physical masses in full QCD are obtained 
from Eq. |l| at p sea = [i ud and [i va \ = [i ud or fx s , 
which are fixed using tt, p and either K or <fi me- 
son masses as input. 

Since our simulations with improved action em- 
ploy a mean-field clover coefficient csw = -P~ 3 / 4 
for the quark action, the leading scaling viola- 
tion is 0(g 2 a). Here g 2 is the renormalized cou- 
pling gj-^((j,) evaluated at a fixed scale fi. 
Higher order contributions consist of terms of 
0(g 4 a log(et)), 0(a 2 ) etc. However, the number 
and range of lattice spacings restrict our choice of 
fitting form to the leading scaling violation term, 
linear in a, which we apply for extrapolations to 
the continuum limit both in quenched and full 
QCD. 

Our main results for the meson spectrum are 
given in Fig. ^ where (j> and K* meson masses 
obtained with A"-input are plotted. The new 

Table 2 

Parameters of quenched QCD simulations. 

16 3 x32 24 3 x48 

(3 a [fm] La [fm] (3 a [fm] La [fm] 

2.187 0.200(2) 3.21(3) 2.416 0.145(2) 3.47(4) 

2.214 0.190(2) 3.05(3) 2.456 0.133(1) 3.19(3) 

2.247 0.181(2) 2.89(3) 2.487 0.128(1) 3.08(3) 

2.281 0.177(2) 2.82(3) 2.528 0.121(1) 2.89(3) 

2.334 0.163(2) 2.61(3) 2.575 0.113(1) 2.71(3) 



quenched results with the improved action nicely 
confirm the finding with the standard action in 
Ref. ||] that the quenched QCD prediction for 
the cj> and K* masses is smaller by about 4% than 
experiment. On the other hand, the two-flavor re- 
sults exhibit a larger slope than those of quenched 
QCD, and the values linearly extrapolated to the 
continuum limit lie within 1% of the experimental 
values. 
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Figure 1. Meson masses in two-flavor (filled sym- 
bols) and quenched (open symbols) QCD. 

We find this result to be a clear evidence of the 
importance of sea quark effects. The observed 
trend that the meson hyperfine splitting is en- 
larged by such effects is consistent with a qualita- 
tive view that the spin-spin coupling in quenched 
QCD is reduced compared to full QCD due to 
a faster running of the quenched QCD coupling 
constant with scale. 

Our baryon mass results are collected in Fig. ^. 
Starting with quenched QCD, our new results 
with the improved action show sizable statisti- 
cal fluctuations. With ten points of data, how- 
ever, one observes a trend that the magnitude 
of scaling violation is smaller than for the stan- 
dard action. Linearly extrapolating to the con- 
tinuum limit, the spectrum for octet and decuplet 
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Figure 2. Corresponding plot for baryon masses. 

baryons is consistent between the two actions. 

A worry with the baryon spectrum is possi- 
ble finite-size effects. As we see from Table § 
however, the spatial size mostly lies in the range 
La rs 2.8-3.5 fm, and only in two cases goes 
down to La w 2.7 — 2.6 fm. Finite-size effects 
are suppressed in quenched QCD due to Z(3) 
symmetry ||. Past studies || also indicate that 
La w 3 fm is sufficient to contain the magnitude 
of finite-size effects for baryons within a few % 
level. We therefore consider that our quenched 
baryon results do no suffer seriously from finite- 
size contaminations. 

Clear conclusions are difficult to draw for two- 
flavor results. Since quenched results are gener- 
ally smaller than experiment, we expect sea quark 
effects to push up the baryon masses. While this 
is a general trend observed in our data, the nu- 
cleon and the A masses turn out higher than ex- 
periment, whereas a better agreement is seen for 
heavier strange baryons. 

Finite-size effects are severer in full QCD with 
dynamical sea quarks since suppression due to 
Z{i) symmetry is no longer effective. It is prob- 
able that the spatial size of 2.5 fm in full QCD, 
which is smaller than the smallest quenched lat- 
tice, is too small to avoid finite-size effects for 



Table 3 

Final results for quark masses in MS scheme at 
H = 2 GeV (in MeV). 







action mud 


m s 
if- input 


^-input 


N f 


= 2 


impr. 3.44±°;£ 


88+ 4 6 


90t 5 n 


Nf 


= 


impr. 4.36±g$ 


not* 


132^ 


N f 


= 


stand. § 4.57(18) 


116(3) 


144(6) 



baryons, especially those made of light up and 
down quarks. Further studies on larger spatial 
sizes are needed to convincingly demonstrate sea 
quark effects in the baryon mass spectrum. 

3. LIGHT QUARK MASSES 

The light quark masses are fundamental pa- 
rameters of Nature which can only be determined 
by hadron spectrum calculations. In Table |^ we 
present our two-flavor results for up-down and 
strange quark masses in the MS scheme at /i = 
2 GeV, and compare them with quenched results. 
These results have been published in Ref . 0] . The 
quenched results, which are consistent between 
the two types of actions, overestimate the quark 
mass by about 25% compared to two-flavor QCD. 
Furthermore, a discrepancy in the strange quark 
mass in quenched QCD depending on the experi- 
mental input is reduced to a level well contained 
in the error of about 10%. Both these features 
demonstrate a significant improvement in the de- 
termination of light quark masses through inclu- 
sion of two flavors of dynamical quarks. 

Our results are obtained from analyses of two 
definitions of quark mass, one derived from the 
axial Ward identity V M A M (x) = 2amf Nl P{x), 
and another m^ WI = (1/k — l/n c )/2a derived 
from the vector Ward identity. For the latter, 
two choices for k c , either the one in full QCD or 
in partially quenched chiral limit, are examined. 
In Fig. [| we illustrate the continuum extrapola- 
tion for the strange quark mass. Lines show a 
linear fit assuming a common value in the contin- 
uum limit, which we take to be our central value. 

A crucial point is the estimation of systematic 
errors, which arise from (i) chiral extrapolation, 
(ii) renormalization factor, and (iii) continuum 
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Figure 3. Continuum extrapolation of strange 
quark mass for improved actions. Results 
from standard action in Ref. [Q arc also shown 
(crosses). 

extrapolation. We estimate the error from (iii) 
by making linear fits separately to each definition 
of quark mass and taking the spread of results. 
The error due to chiral extrapolation (ii) is esti- 
mated by replacing quadratic terms in Eq. [I] by 
terms with power 3/2 which also appear in chi- 
ral perturbation theory and by introducing cubic 
terms in fits of pseudo scalar mesons. 

For the renormalization factors relating the lat- 
tice quark mass to the MS mass in the con- 
tinuum, we employ the perturbative values at 
one loopg given by Z m = 1 + 0.0400 5 |^(/i) 
for VWI quark mass and Za/Z p = 1 + 
0.0308g|^/i) for AWI quark mass. These ex- 
pressions show that the magnitude of one-loop 
corrections is about 10%. One loop correc- 
tions in improvement coefficients have a simi- 
lar size. To estimate further corrections from 
higher order terms, we compare results for 
two definitions of the mean-field improved cou- 
pling constant: flgL(l/a) = (3.648Wi x i - 
2.648PFi x2 )/3/6 - 0.1006 + 0.03149A/, appropri- 

Table 4 



Contributions 


to total 


error in continuum limit . 




stat. 


chiral 


Z-factor 


cont.ext. 




+2.6% 


+1.2% 


+2.3% 


+ 1.7% 




-2.6% 


-2.3% 


-5.0% 


-2.3% 


m s (i^T-input) 


+2.4% 


+1.6% 


+2.2% 


+1.4% 




-2.4% 


-2.2% 


-5.6% 


-2.8% 


m s (</>-input) 


+4.8% 


+1.5% 


+1.7% 


+0.9% 




-4.8% 


-7.6% 


-6.9% 


-1.6% 



ate for the RG-improved gluon action and al- 
ternatively g zjL(l/a) = WixiP/6 + 0.2402 + 
0.031497V/. For the matching scale we use the 
two values fi = 1/a and [i = n/a. 

The systematic errors estimated as above are 
listed in Table [|. We calculate the total error by 
adding them by quadrature. 



4. CONCLUSIONS 

Our first systematic study of Nf = 2 full QCD 
clearly showed the existence of sea quark effects in 
the meson spectrum and quark masses. Further 
progress can be expected in the near future from 
calculations with larger lattices, smaller quark 
masses and, most importantly, inclusion of dy- 
namical strange quark. 
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